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Abstract. - A magnetic particle with atomic spins ordered in an unstable direction is an example 
of a false vacuum that decays via excitation of internal spin waves. Coupled evolution of the 
particle's magnetization (or the vacuum state) and spin waves, considered in the time-dependent 
vacuum frame, leads to a peculiar relaxation that is very fast at the beginning but slows down to a 
nonexponential long tail at the end. The two main scenarios are linear and exponential spin-wave 
instabilities. For the former, the longitudinal and transverse relaxation rates have been obtained 
analytically. Numerical simulations show that the particle's magnetization strongly decreases in 
the middle of reversal and then recovers. 
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' Introduction. — Magnetization reversal in mag- 
netic particles requires transfer of their excessive Zee- 
man and/or anisotropy energy into a heat bath. Mag- 
netic particles are usually considered as single mag- 
netic moments relaxing to external heat reservoirs such 
as phonons, conduction electrons, etc. Very popular 
is a semi-phcnomenological description in terms of the 
Landau-Lifshitz damping and Langevin noise field (see, 
e.g., Refs. [1-3]). The resulting relaxation is exponential 
qr close to it. 

On the other hand, for larger single-domain particles 
that can accommodate internal spin waves, the latter can 
^erve as a heat bath and thus give rise to the particle's 
relaxation and magnetization switching [4-6] . Internal re- 
laxation processes should be much faster than those that 
operate via nonmagnetic degrees of freedom, with impor- 
tant implications for engineering, determining response 
rates of magnetic elements [7] . 

Spin-wave (SW) dynamics in magnetic particles whose 
magnetization changes in time is a new and challenging 
topic. Of course, small deviations from equilibrium, in- 
cluding the case of Suhl instabilities [4,8-11], can still be 
described with the help of a nonlinear SW theory built 
around the ground state. In general, however, one has to 
redefine the spin-wave vacuum by considering SW dynam- 



ics in the frame of the time dependent particle's magne- 
tization [12]. Earlier, a similar approach was suggested 
for the thermodynamics of low-dimensional magnetic sys- 
tems [13] and that of magnetic particles [14]. 

The false- vacuum initial condition can be created by bi- 
asing the system with a magnetic field that leads to the 
disappearance of the metastable energy minimum [5] . An- 
other method is the precessional switching [15] that allows 
to rotate the magnetization by ultrafast magnetic field 
pulses as short as 2 picoseconds with a field amplitude 
as small as 184 kiloamperes per meter. 

Up to date, the coupled motion of the particle's magne- 
tization (the vacuum) and internal spin waves is not yet 
completely understood analytically, and different scenar- 
ios of relaxation are not yet classified. In addition, the re- 
laxation from homogeneous states far from equilibrium has 
not yet been traced numerically over large enough times 
for large enough particles to see all stages of relaxation. 

This Letter presents analytical and numerical work 
showing that the interplay between the global and local 
modes (the magnetization m and spin waves) leads to an 
initially very fast relaxation due to SW instabilities which 
slows down to a long nonexponential tail in its concluding 
stage. Such a kind of relaxation seems to be unknown in 
magnetism. 
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The numericaUy observed magnetization switching is al- 
most complete, and it is typically accompanied by a strong 
reduction of the particle's magnetization m = |m| as the 
spin waves are excited and the system deviates from the 
initial false-vacuum state. Subsequently m recovers to 
slightly less than the initial value 7ti ^ 1 as the excited 
SWs dissipate into the small-amplitude short-wavelength 
excitations and the system approaches the true vacuum 
state with some thermal disordering. This longitudinal re- 
laxation resembles that described by the Landau-Lifshitz- 
Bloch (LLB) equation [2,3]. 

One can identify two main scenarios of the relaxation 
via internal spin waves. These are the linear and exponen- 
tial instabilities, characterized by the corresponding time 
dependences Am(t) = m(t) — 7n(0) at short times. 

The linear instability is a one-magnon process of the de- 
cay of the unstable (false) SW vacuum into the lower-lying 
SW quasicontinuum. In the standard SW theory operat- 
ing near the energy minimum, the one-magnon process 
is of course absent. This mechanism occurs in all cases 
of slightly noncollinear magnetization that include dipole- 
dipole interaction (DDI), surface anisotropy, and random 
anisotropy. The latter will be used below for analytical 
and numerical work as the most transparent case. As the 
magnetization gradually rotates to its equilibrium direc- 
tion, the false SW vacuum becomes closer to the true vac- 
uum and thus less unstable. This leads to a slowing down 
of the relaxation. 

The exponential instability in the case of non-circular 
magnetization precession can be realized in models with 
biaxial anisotropy, uniaxial anisotropy with transverse 
magnetic field. DDI, and many others. Physically, this 
mechanism is a kind of Suhl instability [4,8,9]. However, 
in the global-magnetization frame it finds a much sim- 
pler analytical description [12] that is not limited to small 
deviations from equilibrium. In this case too, numerical 
simulations show a long magnetization tail, after a very 
agile initial relaxation. 

Spin model and relEixation formalism. We con- 
sider a classical spin model (|si| — 1) 

H = Ha - h • ^Sj - i ^ J„Sj • Sj, (1) 

i i'j 

where h = g/i^H is the magnetic field in energy units, Jij 
is the exchange, and T-La is the anisotropy energy. In the 
present work, we consider two models of anisotropy. The 
first is that of random uniaxial anisotropy with the energy 

= -Dr (u^ • s,)' , Dr>0, (2) 

i 

where is a random unit vector. The latter does not 
break the global particle's isotropy, and its role is merely 
to provide amplitudes for SW processes that lead to non- 
conservation of m-h and thus to reversal. The second 
model is that of oriented uniaxial anisotropy, i.e., all spins 



have their easy axes lying in the same reference z direction, 
and thus 

nA = -DY,sl, D>0. (3) 

i 

We define the particle's magnetization as 

™ = ;;^E^»' (4) 

i 

where A^ is the total number of spins. 

The atomic spins obey the Larmor equation 

s, = [s, X n,] , nn, = -dn/ds,. (5) 

We do not include any ad hoc relaxation terms, so that all 
relaxation that we observe is due to the intrinsic nonlin- 
earities in Eq. (1). 

Now, we present the formalism of relaxation via internal 
spin waves using the Hamiltonian (1) with Ti^ being given 
by either (2) or (3). For convenience, we keep both these 
terms in general calculations. The eff'ective field fli in Eq. 
(5) in vector-tensor notations becomes 

= h + 2(D + gO-s, + ^JySj (6) 

3 

where (T))af3 = DSa^z and (g^Q/j = DRUi^aUt^/s, with D 
and Dr being defined in Eqs. (2) and (3), respectively. 

To construct a spin-wave theory around the false vac- 
uum, one has to split s^ into two parts describing respec- 
tively the vacuum mode [the particle's magnetization m 
of Eq. (4)] and the internal spin waves, that is 

s, = m + -0„ ^1/'. = 0. (7) 

i 

Whereas in the standard SW theory m is a constant corre- 
sponding to the ground-state orientation, here it is treated 
as a time-dependent variable. Since the atomic spins are 
subject to the chiral constraint s| = 1, one can use [12] 

m = nyj 1 — ijjf with n • i/ij = 0, where n is a unit vector. 
Although this reduces the number of the ■j/'j components 
in the frame related to the particle's magnetization m to 
only two, the formalism becomes much more cumbersome, 
whereas the final results are not affected. Thus we decided 
not to use the chiral constraint explicitly in our presenta- 
tion. Of course, properly written equations must satisfy 
this constraint that can be used to check them. The equa- 
tion of motion for m following from Eq. (5) has the form 

hm. = [mx (h + 2D • m)] + R, (8) 

where R couples the dynamics of m to spin waves 

^ = -^^{[mx(2g,.,^J] 

i 

+ [V^, X (2D.tA, + 2g,.(m + ^,))]}. (9) 
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In turn, the equation of motion for has the form 



with h <^ J and 



h^.^ = [mx(2g,-m)]+A,|'' +Af\ 



(10) 



where A^^'' is linear in ip^ 



= [ip, X (h + 2D • m + Jom) 



m X 



(11) 



(2) . 2 

and A^ contains terms of order ip and 'ip9- The first term 
in Eq. (10), that causes the spin non-coUinearity, is respon- 
sible for the linear SW instability via direct generation of 
spin waves. We remark in passing that the same effect 
is produced by e.g., surface anisotropy and dipole-dipole 
interactions. The terms containing D in A^^^ can be re- 
sponsible for the exponential SW instability. Since these 
terms are linear in xp^, this instability requires some SWs 

(2) 

in the system to start with. Upon ignoring A^ in Eq. 
(10), one can solve the resulting linear differential equation 
for i/jj and insert the solution into R of Eq. (9). Techni- 
cally, integrating out xp^ can be done most conveniently in 
the particle's magnetization frame with the basis vectors 
(11,61,62), where n = m/m, 6i_Ln, and 62 = [n x 61]. 
This calculation takes into account the time dependence 
m(t) given by Eq. (8) without the perturbation R. Since 
in a large particle spin waves form a quasicontinuum in 
the k-space, integrating over k across the resonance cor- 
responding to the SW instability leads to the part of R 
that describes relaxation of the particle's magnetization in 
Eq. (8). The idea of this calculation is similar to the case 
of spin-lattice relaxation with replacement of phonons by 
internal spin waves, however the implementation is more 
involved since the latter are defined in a moving frame. 

In the case of oriented anisotropy, Eq. (3) or more com- 
plicated forms, the origin of the exponential SW instability 
leading to relaxation is the elliptic precession of spins in 
the spin wave. This complicates the analysis, although 
one can show the existence of SW instabilities [12]. In the 
model with random anisotropy only, D = 0, the depen- 
dence m{t) is a circular precession. Here the calculations 
simplify and lead to the equation of motion for m with 
dissipation due to the linear SW instability that can be 
written in the form 



rh = — [m X h] 

h 



, , \mx [m X hll , ^ m 



m 
(12) 

Here the transverse and longitudinal relaxation rates de- 
pend on the particle's magnetization orientation x = 
(m-h)/(m/i). 

rx(x) = r^o*x(x) = ^^yf'i>x(x) (13) 
r„(^) - r„„$„(.)^^f /fc,„(.) (14) 



$i(x) 
$11 (x) 



1 - X 

6 

(l-x) 



(1 + 2xf + V2{1 -x){2 + x)] (15) 



il + 2xf + V2{1-. 



(16) 



The applicability of our method requires T±, ^ ujh = 
h/h, i.e., the relaxation of the magnetization (as well as 
the rate of SW production) is much slower than the mag- 
netization precession considered as unperturbed in the 
first approximation above. 

We note that Eq. (12) is similar to the LLB equa- 
tion [2,3]. However, it is valid, in general, only for short 
times, when the number of excited spin waves is still small 

(2) 

and m = 1, so that neglecting A] in Eq. (10) is justified. 
Correspondingly, m in Eq. (12) could be replaced by 1. 
More consequent, however, is to introduce the magnetiza- 
tion direction n and write m = nm. Then Eq. (12) can 
be split into the two following equations 



1 r , 1 / N [nx [n X h]] 
n - - [n X h] - T^ixy- L_ ^1, 



where x = n ■ h/h, and 



-T.{x). 



(17) 



(18) 



The small-t behavior mz{t) = 1 — r||( — l)i following 
from Eq. (18) and shown in Fig. 1, agrees well with the 
numerical result, a small discrepancy coming from h/ J 
not being small enough. Eq. (18) does not describe the 
increase of m after switching and its recovery to to = 1, 
that is seen in Figs. 1 and 2. 

A striking feature of our result is that both rj^(x) and 
r||(a;) vanish for m||h (i.e., for x = 1), while they reach 
their maxima at x = —1. Thus, initially fast relaxation 
slows down when the particle approaches equilibrium, see 
Fig. 1. Indeed, Eq. (17) can be rewritten as 



T^{x){l-x^). 



(19) 



In terms of the angular deviation from equilibrium we have 



J/ = 1 - n • h./h = 1 - a:: < 1 



(20) 



upon which Eq. (19) simplifies to y 
the solution 



-3r^o2/^ and yields 



(21) 



l + 3y(0)rxoi 

The long-time asymptote of this solution does not de- 
pend on the initial condition y{Q) and is given by y(t) = 
l/(3r^oi)- The full magnetization vector m = nm in- 
cludes TO that follows from Eq. (18). The latter becomes 
m = -(9/4)r||o2;2 = -(3/4)rio2/' = v/^- This yields 



1 - m(t) = [y(0) - y(0]/4. 



(22) 



One can see that the deviation 1— m(t) remains finite 
and small for t ^ 00, that is the consequence of slowing 
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down the relaxation as m approaches equihbrium. The 
change of m due to its rotation and due to the change in its 
magnitude are comparable with each other. In addition, 
one has to remember that the equation of motion for m 
was obtained for the initial stage of relaxation only and, in 
particular, Eq. (18) does not describe thermalization. Still 
the analytical evidence of the nonexponential relaxation 
via internal spin waves in magnetic particles is sufRcient. 

In the sequel we consider both mechanisms of SW in- 
stability in magnetic particles in more physical terms and 
present the results of our numerical simulations. 

Random anisotropy and linear instability. The 

first particular model we study here is that with random 
anisotropy [16,17] of Eq. (2). In this case, spin waves may 
be generated out of a false vacuum by one-magnon pro- 
cesses. In particular, for m antiparallel to h, spin waves 
energies et in the particle have a negative gap —h. Thus 
a SW can be created if its energy 

Ek = -/i + Jo - Jk (23) 

satisfies the resonance condition 

Ek = 0. (24) 

We recall that k in fine particles is not conserved due to 
the lack of translational invariance. Here Jk is the Fourier 
component of Jij. For a box-shaped particle with free 
boundary conditions (fbc) there are standing spin waves 
with discrete wave vectors [14] 

fca = na=Q,l,...,Na-l, a = x,y,z, (25) 

aNa 

where a is the lattice spacing and N^NyN^ = JV. For 
a particle with simple cubic (sc) lattice one has Jk = 
2J^^ cos(afcQ), so that in the long- wave-length limit 
Jq — Jk = J (afc)^ . For particles small enough, the lowest 
value of Jo — Jk for k ^ [i.e., = 1 in Eq. (25)] may 
exceed h and thus SWs cannot be created. This yields the 
absolute stability criterion for the particle's linear size L 

L < Lmax ^ na\/ J/ h (26) 

in the case of not strongly oblate or prolate particles [c.f. 
Eq. (28)]. Under this condition, there is no relaxation and 
the particle remains in its initial false-vacum state, say, 
antiparallel to h. 

For L > L^jiiix the following possibilities arise. If all 
discrete SW modes in the particle are far from the reso- 
nance with the false vacuum and Eq. (24) is not satis- 
fied, there is no relaxation as before. If there is one mode 
k that exactly or approximately satisfies Eq. (24), this 
mode becomes unstable and its amplitude grows. For in- 
stance, in the sc-lattice particle with = Ny ~ = 5 
{Af = 125) the three degenerate modes with k =(1,0,0), 
(0, 1, 0) and (0, 0, 1) are excited for the magnetic fields in 
the vicinity of /i/J = 0.381966. For the particle initially 



magnetized antiparallel to the magnetic field, the mag- 
netization, obtained numerically, is oscillating because of 
reversible transitions between the false-vacuum state and 
the state with spin waves excited. This phenomenon is 
similar to the probability oscillations between two reso- 
nant states in quantum mechanics. If there are several 
degenerate unstable k- modes and/or the k- modes can be 
converted into other spin-wave modes via nonlinear pro- 
cesses, the time dependence of the magnetization becomes 
more complex and the deviation 1 — increases. Still, 
for small particles it does not lead to a full magnetization 
reversal since the spin-wave modes are discrete and it is 
difficult to fulfill the energy conservation laws in nonlinear 
SW processes. 

To illustrate the linear instability in larger particles 
with quasi-continuous SW spectrum, we have solved Eq. 
(5) for magnetization switching in the random-anisotropy 
model. The particle is a box with fbc, consisting of 
28 X 31 X 34 = 29512 atomic spins. In Fig. la we plot the 
results for an initial magnetization antiparallel to the field. 
One can see a nearly full reversal that goes linearly at small 
times. For this initial condition, the transverse magneti- 
zation component m± remains small. In the middle of 
the switching, excited spin waves almost destroy the mag- 
netization m. Then these long-wavelength SWs of high 
amplitude convert via nonlinear processes into all possible 
spin-wave modes, and the system thermalizes. Since the 
thermodynamics of classical spins is mainly determined by 
short-wavelength modes of high energy, the energy conser- 
vation requires that the amplitudes of these SWs be small. 
This explains the almost full recovery of m after reversal. 
The asymptotic disordering 1 — m > (see Fig. 1) is due 
to the final temperature T corresponding to the released 
Zeeman energy. The numerical results in Fig. lb (initial 
magnetization perpendicular to the field) show that the 
magnetization relaxes mainly by rotation towards equilib- 
rium, although m decreases from saturation in the middle 
of the process and then recovers. 

Oriented anisotropy and exponential instability. 

— The second model of relaxation via internal spin waves 
we consider is the model with oriented uniaxial anisotropy 
(3). In this case, the instability turns out to be exponen- 
tial. The latter occurs for a general form of the single-spin 
part of Eq. (1) that leads to non-circular magnetization 
precession, see Ref. [12]. For illustration, we use Eq. (3) 
in (1) with h ^he^ and we choose the initial state s^ = e^- 
Linearization around the latter yields the spin-wave spec- 
trum 

£k = Vih + Jo - Jk) ih-2D + Jo- Jk). (27) 

For h > 2D the state s^ = e^^ is the energy minimum and 
£k is real. In the interval < h < 2D the state Sj = e^, is 
the energy saddle point, and m will rotate away from this 
state, similarly to the case studied in Ref. [5]. SW modes 
in the interval < Jo — Jk < 2D — h are unstable since Ek 
is imaginary [12]. For h < the state s,; = e^, is the energy 
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Fig. 1: Magnetization switching via linear spin- wave instability in a particle with uniaxial random anisotropy. a) Particle 
prepared with all spins antiparallel to the magnetic field. The small-t asymptote m^{t) = m{t) = 1 — T\\[—l)t of Eq. (14) is 
shown by the dashed line. For this initial condition, switching occurs predominantly via changing the magnetization length m. 
b) Particle prepared with all spins perpendicular to the magnetic field. 



maximum, so that m performs small precessions around 
Bx in the absence of SW processes. Indeed, according to 
Eq. (27), the mode k = is stable, so that the particle's 
magnetization cannot rotate away from this state. On the 
other hand, the modes with k ^ in the interval —h < 
Jo ^ J'k < 2Z3 — h are unstable. Thus this field range 
the only way of reversal is via excitation of internal spin 
waves with k 7^ that strongly reduce the magnetization 
magnitude m = |m|. 

In magnetic particles, SW modes are discrete, and the 
instability requires existence of SW modes inside the insta- 
bility interval. The particle is stable if the smallest value of 
Jo — Jk with k 7^ (i.e., Ua = 1) exceeds the right bound- 
ary of the instability interval, i.e., 2D — h. This leads to 
the stability criterion 



imax = aiVmax < TTa^ ^ (28) 

for the particle's size imax = niax(L2:, Ly, Lz)- Eq. (28) is 
similar to the single-domain criterion since its right-hand 
side is the domain-wall width, if /i = 0. 

Efficient magnetization reversal via the SW instability 
requires a quasicontinuous SW spectrum, i.e., many SW 
modes in the instability interval. For this, a strong in- 
equality opposite to Eq. (28) should be fulfilled, i.e., the 
particle's size and/or the anisotropy D should be large 
enough. The numerical results for the particle of the same 
shape as above are shown in Fig. 2. The applied magnetic 
field has been chosen stronger than the anisotropy field, to 
make the equilibrium magnetization parallel to the field, 
for the sake of comparison with Fig. 1. For the initial 
magnetization antiparallel to the field (Fig. 2a), the first 
stage of the relaxation is an exponential increase of 1 — m 
due to the SW instability until the first dip, followed by 
an incomplete recovery of m and its further decrease. The 



subsequent development is similar to that in Fig. la. For 
the initial magnetization perpendicular to the field (Fig. 
2b) the behavior is more complex as the components of m 
perform fast oscillations because of the elliptical preces- 
sion. Apart from the latter, the results at not too small 
times are qualitatively similar to those in Fig. lb, the main 
feature being a decrease and a recovery of m. The latter 
was observed in our simulations for other orientations of 
the field with respect to the anisotropy axis. 

Some orders of magnitude. — Let us now estimate 
the magnetization reversal rates due to the internal SW 
excitation. First, since Eq. (1) is written in terms of the 
normalized spins = 1, conversion to actual spins S 
requires replacements J => JS and D =^ DS, while h 
remains the same. For cobalt S = 3/2, and the above 
modification is inessential for estimations. Comparison 
of Figs. 1 and 2 shows that the exponential instability 
leads to a faster relaxation than the linear instability. In- 
deed, for all materials D ^ J is satisfied {D ~ 0.0024J 
for Co), and the instability increment in Eq. (27) is of 
order D/fi, while Eq. (13) contains additional small fac- 
tors D/J and y^h/J. However, the difference of the over- 
all relaxation rates is not that big since at larger times 
the relaxation rate should be of the order of in both 
cases. Since we do not have any analytical result for the 
exponential scenario at large times, we use Eq. (14) for 
estimations. With J/ks ~ 10^ K for Co (assuming amor- 
phization leading to random anisotropy), one obtains the 
rate F w 10^ s~^ for h ~ O.OOIJ that corresponds to 
H w 0.07 T in real units. For this field, the precession 
frequency is luh = {gf^s/^)-^ ~ 10^° i-'^-j precession 
is fast with respect to relaxation. Preparing the particle's 
magnetization in a desired direction requires the applica- 
tion of pulsed fields [12, 15] that are at least by an order 
of magnitude larger than H, during a time that is at least 
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Fig. 2: Magnetization switching via exponential spin- wave instability in a particle with oriented uniaxial anisotropy and trans- 
verse field, a) Particle prepared with all spins antiparallel to the magnetic field. Again the longitudinal relaxation is the main 
mechanism of switching, b) Particle prepared with all spins perpendicular to the magnetic field. 



by an order of magnitude shorter than 1/ujh- For the ex- 
ponential scenario, one can expect a similar result for F 
without the factor ^Jh/J, i.e., F « 3 x 10^ s~^. The pre- 
cession frequency due to the anisotropy field oja in Co 
is comparable with ujh quoted above. These estimations 
give an idea of the relaxation rates that remain in the 
same ballpark for different kinds of the anisotropy. For 
instance, the linear instability can also be driven by the 
surface anisotropy that is much stronger than the volume 
anisotropy but acts on the surface spins only. Numerical 
simulations with surface anisotropy Ds <C J give the re- 
sults similar to those for the random anisotropy. If the 
surface effects are so strong that another type of ordering 
is realized at the surface, the order in the core should be 
strongly distorted as well. This will presumably lead to a 
very strong internal relaxation, although simulations are 
needed to clarify the situation. 

Our numerical simulations for a multitude of different 
models and initial states of the particle use exaggerated 
values of the anisotropy and field. Using realistic smaller 
anisotropy and field values leads to the necessity to work 
with much larger sizes, according to Eqs. (26) and (28). In 
addition, the relaxation becomes too slow, so that much 
longer time interval is needed for computations. Both fac- 
tors together make this impossible. Nevertheless, for par- 
ticle's sizes M > 10** the results for realistic parameters 
can be obtained from those for the exaggerated parame- 
ters by rescaling. 

Conclusion. — We have shown that exponential and 
linear spin-wave instabilities in large enough magnetic par- 
ticles put into a false- vacuum state (i.e., with the mag- 
netization deviating from the equilibrium direction) can 
lead to complete magnetization reversal and thermaliza- 
tion. Whereas the former is similar to the Suhl instabil- 
ity [8,9], the latter seems to have no analogue within the 
standard SW theory. It should be stressed that both kinds 
of instabilities are described within the linear approxima- 



tion in spin-wave amplitudes, unlike the description within 
the standard (nonlinear) SW theory. In both cases, the 
relaxation is fast at the beginning but then slows down. 
These results should be useful for investigating the dy- 
namics of magnetoelectronic elements and, more general, 
in the physics of unstable macroscopic states. 

* * * 

D. G. thanks E. M. Chudnovsky for useful discussions. 
D. G. is a Cottrell Scholar of Research Corporation. 
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